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ABSTRACT 

We calculate the creation of electron-positron pairs in Gamma- Ray Bursts (GRBs) 
resulting from the collision between scattered and outward moving gamma-ray pho- 
tons. The number of pairs exceeds the number of ambient medium electrons encoun- 
tered by the GRB ejecta up to ^ 10^^ cm from the center of explosion. The shock 
resulting from the interaction of the ejecta with the pair-wind may brighten the af- 
terglow synchrotron emission during the first few minutes. Even without this effect, 
the peak intensity of the optical afterglow increases with the density of the surround- 
ing medium. Therefore, observations of the optical flux at early times constrain the 
density of the circumburst medium. If the electron and magnetic field energies be- 
hind the forward shock sweeping-up the pair-wind and the circumburst medium are 
as inferred from fits to the broadband afterglow emission at 0.5 — 100 days, then the 
current upper limits on the optical counterpart emission, set by the ROTSE and LO- 
TIS experiments, indicate that the circumburst medium within 0.01 pc is less dense 
than 100 cm"'^ or, if a wind, corresponds to a progenitor mass-loss to wind speed ratio 
below lO-6Af0/yr/(lO3km/s). 
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1 INTRODUCTION 

The 7-photons produced in Gamma-Ray Busts (GRBs) are 
scattered by electrons in the medium surrounding the ex- 
plosion, their interaction with the outward moving photons 
creates an electron and positron pair if the photons energy 
in their center of mass frame exceeds 2meC?. These pairs 
scatter more photons, leading to an exponential growth of 
their density. The pairs are accelerated by scattering the 
gamma-ray burst photons. 

The process of pair-creation in GRBs has been previ- 
ously considered by Thompson & Madau (2000), Meszaros, 
Ramirez- Ruiz & Rees (2001), and Beloborodov (2002). The 
resulting pair-wind moves at a relativistic velocity and the 
pairs dominate the number density of charged particles up 
to a distance 10^® (iJ^/10^^)^'^^ cm from the center of ex- 
plosion, where is the isotropic-equivalent GRB energy 
release. The purpose of this paper is to consider the effect 
of pairs created by GRB pulse on early afterglow emission. 

In §2 we consider all of previously discussed physical 
processes for electron-positron pair creation by GRB pho- 
tons. We include in our calculations the overtaking of pair 
plasma by the material ejected in the explosion, which can 



substantially reduce the number of pairs produced. In §3 we 
consider the effect of pairs on early afterglow emission, and 
calculate for various external densities and ejecta Lorentz 
factor the optical synchrotron emission produced by the 
shock that heats up the pair-wind when it collides with the 
GRB ejecta. The non-detection of early optical afterglow 
emission brighter than magnitude _R ~ 13 by the ROTSE 
and LOTIS experiments can be used to provide an upper 
limit on the density of the circumburst medium. 



2 DYNAMICS OF THE 

ELECTRON-POSITRON WIND 

The enrichment is calculated from the energy spectrum 
of the incident and scattered photons. We approximate the 
GRB energy spectrum as a power-law F(e) = Fp{e/ep)~" 
for e > Ep. For most of the 5,500 spectra of the 156 bright 
BATSE bursts analyzed by Preece et al. (2000), the spectral 
slope a above the Sp break is between 0.5 and 2.5, with the 
peak of the distribution at 1.3, while Sp is between 70 keV 
and 800 keV, its distribution peaking at 250 keV. Assum- 
ing a burst redshift z = 1, the corrected break energy Sp, 
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measured in rneC^ units is between 0.3 and 3, with a peak 
at 1. 

The flux of scattered photons at some distance x from 
the leading edge of the gamma-ray pulse, is calculated nu- 
merically by integrating over the scatterings that occurred 
from X = to the current x and over the primary photon 
spectrum and deflection angle, taking into account the pho- 
ton attenuation during its propagation. 



2.1 Planar Geometry 

In planar geometry, the laboratory frame number density 
n± of the created pairs and their streaming velocity /3c are 

given by 



(1) 



{ [2n±me + np(me + m^)] 7/3(1 - 0)c''} = ,(2) 



dx 



dt 



respectively, where t is time measured in the laboratory 
frame, Up is the laboratory frame number density of the 
swept-up external medium electrons, rUp is the proton mass, 
and 7 is the pair-wind Lorentz factor. The right-hand side of 
equation (1) represents the rate of pair creation from scat- 
tered and primary photons, while the two terms in the rhs of 
equation (2) represent the rate of momentum deposition by 
pair formation and by photon scattering, respectively. Con- 
servation of the number of external electrons implies that 



d_ 
dx 



[np(l-/3)] = 0. 



(3) 



i.e. Up = 7i/(l — 0), n being the density of the unperturbed 
circumburst medium. 

The most general form of the rhs's of equations (1) and 
(2) are given by Beloborodov (2002) (see his equations [10] 
and [16]). For analytical purposes we simplify them by mak- 
ing two approximations. First, we ignore the angular depen- 
dence of the photon scattering probability and consider that 
an electron (or positron) at x deflects an incident photon 
only by the average deflection angle cos~^/3(a;). This ap- 
proximation reduces the integral of the scattered flux from 
one over the wind volume to one over a surface from which 
scattered photons arrive at the same time at x (where 
creation takes place). Secondly, we consider that the e* pairs 
are produced mostly by photons near the peak of the spec- 
trum, i.e. in the integral over the spectrum of primary pho- 
tons wo approximate by unity the term 1 +0/7, where e 
is the photon energy in rrieC^ units. This approximation is 
justifled by that photons with e > 7 are scattered in the 
Klein-Nishina regime, thus the scattering probability is re- 
duced. Taking into account that the GRB photon spectrum 
above meC^ is typically steeper than £~^, the latter approx- 
imation is more suitable at locations x where the wind is 
relativistic (7 » 1). 

With the above approximations, equation (1) for the 
wind number density becomes 



dx 



[n(l-/3)] 



(a) 



A2 



X 

jdx' 



h{x) 



■Pix') 



l + Pix') 



f /(£)/(£*) 



/(-) = ( r 



(4) 



where n = 2 n± -\- Up is the total lepton density, and <j>{oi) 
(7/12)(l-|-a)-'5/3 (Svensson 1987), 



Fp(JT ' 

ctt being the Thomson scattering cross-section, and 

2(1 + /3) 



£fc„ ~7(l-F/3) 



(5) 



(6) 



£(1-/3) 

are the photon threshold energy for formation and the 
Klein-Nishina "cut-ofF' photon energy, respectively. Note 
that for £p < 2 (i.e. the peak of the GRB vFu spectrum 
is below ~ 1 MeV), £* is larger than Ep. 

With the same approximations, equation (2) for the 
wind momentum density becomes: 



d_ 

dx 



n+^np 7/3(1-/3) 

rUe 



n{l - fi) 
A(l+/3) 



/ 



def{e) 



+ 4^/d.'n(.')i^ / f /(.)/(..) P., 



where 

p± = ^ + ^^^i+^ 3' = 3{x') 



(7) 



(8) 



is the total radial momentum of the scattered and primary 
photons that form a pair. Equations (4)- (7) arc similar 
to equations (42)- (47) of Beloborodov (2002). 

For a power-law photon spectrum, equation (4) leads to 



^Ki-/?)] = ^^ 



[(1+/3')7T' 



'.a+2 



The integrand in the above equation decreases rapidly with 
7, most of the contribution to the integral coming from 
7^2. Physically, this means that pairs are created by pho- 
tons scattered by mildly or non-relativistic pairs; the 
threshold for photons scattered by relativistic is much 
greater because these photons are moving within an angle 
7^^ relative to the primary photons. 

As shown in Beloborodov (2002) the momentum de- 
position in the e^ wind through pair creation is at most 
comparable to that transferred by photon scattering. For 
simplicity, we retain only the scattering term in equation 
(7), reducing it to 



d_ 

dx 



{hme + npmp)j3 
7(1+/?) 



me£p n 
(a-l)A (l-h/3)V 



(10) 



2.1.1 Non-Relativistic Wind 

In the /3 <gC 1 limit, equations (9) and (10) lead to an expo- 
nential growth of the pair density and velocity: 



h{x) = ne"/^± , p{x) = 



X± x/x + 

e ' * 



(a — 1) nip A 



(11) 
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for X > x±, where 



x± = 



A a = 1.5 



-3/2 



0(a) eg ^' 

is the length-scale for pair loading. For a GRB of isotropic- 
equivalent output energy E~f, equation (5) gives for the ac- 
celeration length-scale at radius r 



A 



1 V 0(a) OT e^-'-Et, 



0.1 



rn 1/2 
-C'7,53£p 



, (13) 



where A = cT/{l + z) is the geometrical thickness of the 
GRB front, T being the observed burst duration, z the 
GRB's redshift, and the usual scaling X„ = X/10" was 
used. 

Prom equation (11) we see that the wind acceleration 

longth-scalc Xacc, defined by 7(x-acc) = 2 is a factor of a 
few larger than x±. A simple estimate of Xacc can be ob- 
tained by equating the momentum deposited by photon scat- 
tering hacc{cTTXacc)(Fpepme(? /(?) (mec)(a;acc/A) 

(considering that all photons have the same energy meC^Ep) 
with the momentum [naccTrie + np^accmp)c ~ nrUpC of the 
medium (assuming that the wind mass density is dominated 
by the protons, i.e. n{xacc) <S np{xacc){iTip/me), and that 

^p(^acc) ^ ^)- 

£p = l 

-* Xacc » 5X± . (14) 



''acc ^xacc/x+ 



A 



rueep 



Then equation (13) shows that the wind is accelerated to a 
relativistic speed for radii smaller than 



rend " W''e'J'E!^% cm . 



2.1.2 Relativistic Wind 



(15) 



For X » Xacc we can solve the pair density and momentum 

equations together to obtain h{x) & 'y{x). We consider a 
power-law solution of the form 

ni,^=Aj-^y\ ^ix) = ^(^y\ (16) 

\>^acc^ \.*^accy 

Substituting this into equations (9) and (10), and assum- 
ing that the mass density is dominated by protons i.e., the 
rmmbor of pairs per proton is much less than mp/nie, we 
obtain 



(ai — 202)^2' 
a2nmpy°'^~^ 

Xacc 



„ ai_2o2-l „ 2 r /»l-2(" + l)°2 

y =27] y 



(17) 



nicEp 



^1 \ ai-2a2 



(18) 



where y = x/xacc, ij = Xacc/x± is a constant independent 
of X, and Xmin is determined by the condition that the pair 

creation optical depth for a photon scattered at Xmin, of en- 
ergy ~ £p, is unity between Xmin and x; Xmin 's determined 
by the following equation 



(^X ^mzn) 



£t 



2(1 + '^) -^.(19) 



By substituting (16) in the above equation we find Xj, 
explicitly 

2a 



2«A 



27 



2°'XA; 



p j ^acc 
■^9 \'^rmn 



1, (20) 



ivy) 



l/(2ao2) 



or 

•^acc ^ 

Xmin ~ yl2''"^[0(Q:)2«+l]V(4aa2) ■ 

It follows from equation (18) that 
3a2 = oi -h 1 



(21) 



(22) 



Furthermore, there are two cases to be considered for the 
integral in equation (17) - one of which is Oi — 2(q -|- 1)02 > 
—1. The integral in this case is dominated by the upper 
limit, and by equating the exponents of y on the two sides 
of the equation we find aa^ = 1. This together with equation 
(22) implies that a < 1/2, which is not a case of interest for 
GRBs since the high energy spectral index for GRBs has 
Q > 1. The other possibility is that ai — 2{a + l)a2 < — 1, 
and in this case pair production is dominated by photons 
scattered at the smallest x i.e. at Xmin. Making use of the 
expression for Xmin (eq. 21) in equation (17), and setting 
the exponent of y to zero results in 



Ol — 202 



oi — 2a2 + 1 
2Qa2 



(23) 



Solving for oi and 02 from equations (22) and (23) we find 



01 = 2-|-3/(2a), 02 = 1 -h l/(2a) 

In other words the solution for n and 7 are given by 



(24) 



h{x) = Ai 



2+^ 



A2 



\Xacc / 



, (25) 



with Ai ~ nexp(?7), the density at Xacc, and ^42 ~ 4. A 
somewhat more accurate expression for Ai and A2 can be 
obtained by combining equations (17), (18) & (24). The re- 
sult in Beloborodov (2002) is a special case (a = 1) of the 
solution presented here. 



2.2 Spherical Expansion and Saturation of Wind 
Acceleration 

Equations (4) and (7) were derived under the assumptions of 
planar geometry and stationary solution in the x-coordinate. 
Due to the spherical expansion, the rate of momentum de- 
position decreases as r^^, while the rate of pair creation 
decreases as r~*. The effect of spherical expansion is easily 
accounted for in a numerical treatment by decreasing the 
photon flux: F{r) = F{r-y){r /r~/)'^ , where r-y is the radius 
at which the GRB photons were emitted, and by adjusting 
similarly at each timestep the laboratory frame density of 
the created pairs n±{x) and accelerated external electrons 
np{x). 

The Lorentz factor jsph and pair density fisph at which 
the pair creation process saturates duo to spherical expan- 
sion {photon dilution) can be calculated approximately as 
the value of these variables at the coordinate Xsph corre- 
sponding to a twofold increase of radius r at which the exter- 
nal medium entered the GRB front: r = J^""" dx' {1- P')'^ . 
Using equation (25), one obtains 



Xsph '■ 



and 



{3 + a-^)r 



8Xa 



a/(l-|-3a) 



(26) 
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^sph ■ 



(3 + a- 



2 



(3 + Q~^)r 



4q + 3 
2(l + 3a) 



1 + 2q 
2{l + 3c«) 



where 



1.5 X 10-^r?,i 



(27) 



(28) 



(29) 



At X > Xsph, the wind density decreases as n oa R ^^'^ cx 
X is the front radius corresponding to 



^ Isph 



2(l + 3a) p 
3+8a 



2.« X iO--^i:\.r,:!(:! + a ^) 



T(l + ^)-ir,ep °{<?i(a)2i-«}-i/2 



l + 2e 
3+8c 



Thus, if the GRB emission is produced in internal shocks 
occurring at ~ 2 V'^cSt - 6 x lO" Tl5t-i cm, where 5t ~ 
0.1 s is the GRB variabihty timescale, the lepton density 
and Lorentz factor of the medium in the wake of the GRB 
front, at r ^ 10^*cm and a = 1.5, are 



ne ~ 3 X 10^ 



7e - 170 



ri4 
T25t- 



(31) 



Tlsph — -^1 



2.8 X 10^(3 + Q-')i;^,53VWa)2^ 
T(l + 0)-iriBr7£-" 



Isph '■ 



2 



2.8 X 10^(3 + a-^)S^,53 V0(a)2^ 



r(l + 0)-iri5r7ep' 



3+4o 
2 + 6a 



1 + 2q 

2 + 6q 



(£T,53£p/rri5)^'''^+'*"^ If a^sph > A the traihng edge of the 
GRB front is encountered before the wind acceleration and 
pair loading is saturated. This happens for electrons entering 
the GRB front at radii larger than 



^ 10^^(6.6 X 10^)^ [r/(i + z)Y^ 
[2.8 X 106(3 + a-i)]^ 
In this case, substituting x 



-E-y,53£p 



l + 2a 
2+5q 



cm. (34) 



x~ , where 7? ■ 
coordinate x. 

The Lorentz factor of pair-wind is also limited by the 
non-zero transverse momentum of gamma-ray photons: the 
angle 6 of a GRB photon w.r.t. the radial direction is 
9 = r~^(r-y/r), where T is the Lorentz factor of the rel- 
ativistic ejecta which produced the GRB photons at r^. 
Consequently, when the pair-wind reaches the Lorentz fac- 
tor 79 = r{r/r-y) the photon field becomes isotropic in the 
frame of the moving scattoror and there is no further net 
transfer of momentum from photons to the wind. As shown 
by Beloborodov (2002), this effect yields a multiplying fac- 
tor [1 — (7/73)*] in the Psc term in equation (2), where 
73 = r(r/r^). From equation (25), the acceleration satu- 
rates at coordinate xe/xacc — (2rr /A2r^)^°'^^^''^^°'\ where 
the wind density is he ~ yli(2rr/A2r^)'^+*"^/(^+2">. 

Using the above equations for r{x), Xsph, and xe, it can 
be shown that the angular dispersion of GRB photons is 
more important than the photon dilution effect (i.e. xe < 
Xsph) at radius r less than 



nA — n exp(?7) 



A in equation (25), yields 

^ 3+4o 



6.6 X 10^ £^,53^0(0)21 



and 



7A 



A2 
2 



6.6 X 10'^E^,r,.iy/(t){a)2^ 
r-knep" 



(35) 



(36) 



for the maximal wind Lorentz factor and pair density, 

thus the pair-enrichment factor is /xa = ftA/(2n7A) oc 

(e^B^,53/rf5)'/'"- 

The main results of this section is that the pair cre- 



ation and acceleration up to r,. 



10^ 



cm is limited by 



the decrease in photon density due to spherical expansion. 
The pair creation for r** < r < rend ~ 10^® cm contin- 
ues throughout the gamma^ray pulse and the results in this 
regime are the same for spherical and plane parallel geome- 
tries. The medium beyond rend is not affected by the passage 
of the GRB front. 

Most of the kinetic energy of the pair-wind is acquired 
at r ^ rrei\ Vrei - 2.6 X 10i''(S,,.53v/#^)2i^ep/77)'/^cm, 
is the radius whore the pair-wind Lorentz factor drop be- 
low 2. Since the pair-enrichment factor jj, does not exceed 
(30)mp/me, the wind kinetic energy is approximately that of 
the entrained circumburst medium protons (or ions) . Using 
equation (36), the wind kinetic energy is found to be 



2 

r 7A 



dr ^ 10''^n£p"''^iS;^^g^3 erg 



(37) 



and therefore the fraction of energy in gamma-ray pulse 
converted to the creation of a pair- wind is 10"'^no. The 



As shown in equation (3), the acceleration of the external 
electrons compresses them to a density np,e = 2 7|n, there- 
fore the pair-enrichment factor /i = h/up (i.e. the num- 
ber of e" and created for each external electron) is 
He ~ 500 riA/{V26t-i). 

The medium encountered at r > r, is pair-loaded and 
accelerated up to the density and Lorentz factor given by 
equations (27) and (28). By substituting Xacc{r) from equa- 
tions (14) and (13), we obtain 



overall wind Lorentz factor 7^ 



3, where 



iriw ~ (47r/3) rg„^nmp is the wind mass. 



(32) 



.(33) 



Note that A\ = hsph{x = Xacc) — nexp(r]), and A2 ~ 
4. The pair-enrichment factor Hsph = "■sph/(2n7sph) oc 



2.3 Numerical Results 

There are two factors associated with the wind spherical 
expansion which have been ignored in the analytical calcu- 
lations presented in section §2: the spherical dilution of the 
incident GRB photon flux and the interactions within the 
pair-wind. The former has been considered in §2.2; here we 
note that another way in which it leads to a non-stationary 
h{x) and 7(2;) is that the GRB front radius may increase 
significantly while a shell of external matter, undergoing ac- 
celeration and pair-enrichment, is within the front. Ignoring 
this effect is appropriate when the upper limit to the ra- 
dius increase during the front-crossing, 27aA, is less than 
the front radius, r, which, after using equation (36), leads 
to r > 2 X 1015^016 cm for ^ ^ ^ iqB3 gj.gs, g^, = 1. 

Collisions within the pair-wind occur because the external 
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10^ 
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\ 



10' 
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(cm) 



Figure 1. Wind acceleration and pair-loading for a GRB pulse 



of isotropic-cquivalent energy = 10 erg, intrinsic duration 

c^, post-peak slope 
is assumed to arise 

at r-y 10^* cm from electa moving at F = 100. Calculations 



10 s and peak of uFi, spectrum Sp 
1.5 (F^ oc 1^-°'). The GRB emission 
y 10^* cm from ejecta moving at F : 



take into account the interactions within the wind, between the 
wind and the GRB ejecta, and the effect of spherical expansion on 
the incident flux and wind density. The abscissa gives the radial 
coordinate at which the ejecta-wind shell interactions take place. 



medium entering the GRB front at smaller radius is accel- 
erated to a higher Lorcntz factor and may overtake a slower 
shell of external matter entering the GRB front at a larger 
radius. A third factor is the presence of the ejecta within the 
GRB pulse, which sweep-up the pair-wind and may drasti- 
cally reduce the acceleration through photon-scatterings if 
the ejecta Lorentz factor is significantly larger than that 
of the wind. Thus this factor is of importance only when 
the lag 5 between ejecta and the leading edge of the GRB 
front, r/(2r^), is less than the pulse width. A, i.e. at radii 
r < 6 X lO^^r^Ti cm, and only if T > j{S). 

The effects discussed above are taken into account in 
the results shown in Figure 1. The external medium is dis- 
cretized in shells of thickness much smaller than their radius, 
which entering the GRB front at a; = and are accelerated, 
compressed, and pair-loaded as described in §2. To include 
the spherical expansion, the GRB flux and the shell density 
are adjusted after each time-step. The progressive merging 
of shells (interactions within the wind) is accounted for by 
replacing a pair of colliding shells with one whose Lorentz 
factor and density are calculated from energy conservation 
and the jump conditions at shocks, respectively. The loca- 
tion of the ejecta is tracked, allowing for the accumulation of 
shocked wind shells and compression of the shocked ejecta. 

The termination of the wind acceleration while it is still 
within the GRB front, due to the collision with the ejecta, 
is illustrated in Figure 1 by the rising part of the wind den- 
sity (n±) and Lorentz factor (7±) at r ~ 10^* cm. As the 
ejecta lag more behind the front's leading edge, the external 
medium is more accelerated and pair-loaded. The interac- 
tions within the wind and the decrease of the incident flux 
with radius lead to a flattening of the wind Lorentz factor at 
r = 2 — 4 X lO'^'* cm, followed by a decrease. The ejecta drive 



the forward shock into a relativistie wind up to r^ei ^ 10^^ 
cm, thus the dissipation efficiency of this shock is reduced 
by the existence of the pair-wind up to an observer time 
Uei ~ (1 + z){rr-ei/cr^) ^ lOFj ^ s. The pair production 



becomes negligible at 
observer time t„a 



~ 10 cm, corresponding to an 
lOtrei ^ 100r2 ^ s. As shown in Figure 



1, the pair enrichment factor decreases roughly as /u oc r 



for 10^^ < r < 10' 
ating leptons at r 



cm. This means that most of the radi- 
= rpair are the pairs formed at r < rpair 



and not the electrons originally existing in the cireumburst 
medium. Then we expect the lepton enrichment to have an 
effect on the afterglow emission until later than tpair- 



3 EARLY AFTERGLOW EMISSION 

A shell of cireumburst medium that enters the GRB front 

at radius Tin leaves the front at radius rout ~ Tin + A/(l — 
f3±), where /3± is the shell velocity reached when it exits the 
GRB front, the corresponding Lorentz factor y± = 'f{x = 
A) being given by equations (31), (33) or (36). The GRB 
ejecta collide with the pair-wind at Tc [1 + {j±/T)'^]rout 
(assuming that 7± <iC F) if the GRB ejecta are outside the 
front, i.e. rc > 2 F^ A. If the GRB ejecta interact with the 
pair- wind while the layer is still within the GRB front, its 
pair-loading and Lorentz factor are lower. 

In the frame of the unshoeked pair-wind, the GRB 
ejecta move at the Lorentz factor jr ~ F7±(l — /3±/3e), where 
/9e is the ejecta velocity in the laboratory frame. Denoting 
by Se the fractional energy of the shocked wind in e~ and 
e"*", the typical lepton random Lorentz factor is: 



p — 1 \ nrrie ) 



(38) 



where p is the exponent of the power-law energy distribution 

dn/''fc oc 7~'' at 7e > 7p acquired by the leptons. 

At Tc the wind density is smaller by a factor {rc/rout) 
than that at rout- The shock compression at Vc increases 
the comoving density of the wind by a factor 47,. -|- 3. If 
the magnetic field stores a fraction es of the shocked wind 
internal energy, then the magnetic field is: 

= 87res(47r.-|-3)(7T.-l)(n±me-|-np,±mp)(ro«t/rc)^ ,(39) 

where n± = n{x = A) and np,± s np{x = A). 

Figure 2 shows the evolution of the characteristic syn- 
chrotron frequency Vp = (3e/16mec)7p-BF corresponding to 
a typical electron energy in the shocked wind (eq.[38]). As 
discussed in §2.3, the pair-wind affects the afterglow emis- 
sion by reducing the energy per lepton in two ways: until 
trei ~ 10 Fj^'^ s, it lowers the dissipated energy and, un- 
til some time later than tpair ^ 100 F^'^ s, it increases the 
number of particles among which this energy is shared. Both 
these effects lower the typical lepton Lorentz factor fp, lead- 
ing to a softening of the afterglow emission. As the GRB 
front radius increases, the wind is less accelerated and pair- 
enriched and the afterglow spectrum hardens, its peak fre- 
quency Vp asymptotically reaching a constant value. Note in 
Figure 2 that the softening of the spectrum depends on F, 
for higher initial Lorentz factor, the gap between the leading 
fronts of the GRB pulse and ejecta is smaller, thus the pair- 
wind is less accelerated and enriched before it is swept-up. 

When the mass of the swept-up cireumburst medium 
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Figure 2. Evolution of the peak frequency of the synchrotron 
spectrum i/p for three sets of parameters (n, F) (external particle 
density and fireball initial Lorentz factor). The burst is placed at 
redshift 2 = 1 (redshift), other parameters being as for Figure 
L The kinetic energy of the ejecta f?fe was talcen equal to 
(i.e. a 50% efBciency of the GRB emission was assumed). The 
power-law energy distribution of the shock-accelerated leptons, 
dn/dje oc Je"' , has a total energy equal to a fraction 6e = 0.1 of 
the post-shock energy and an index p = 2.5. The magnetic field 
energy is a fraction sb = 10^* of the internal energy. The photon 
arrival time is calculated for the photons moving along the line 
of sight toward the center of explosion. Note that the pair-wind 
softens the afterglow spectrum up to 10 seconds. 

reaches a fraction ~ of the ejecta mass, the fireball 
begins to decelerate. This happens at 

-'=(3;^)"'~""(Sl)"'™' 

where Ek is the kinetic energy of the ejecta after the 
GRB phase, corresponding to an observer time tdec = 
(1 + z)ri/{r^c) ~ 700(Ek,53/noy^^r-^^'' s. At t > Uec, 
the fireball Lorentz factor is F = [Ek/m{r)]^'''^ (for adia- 
batic dynamics), where m(r) is the mass of the energized 
circumburst medium. Due to the decreasing F, the after- 
glow spectrum softens again, the peak frequency evolving 
as 

~ 3 X W'^El%el_,e]^%^^/^ Hz . (41) 

To obtain the afterglow optical light-curves, we in- 
tegrate the synchrotron emission function over the elec- 
tron distribution, taking into account effect of the spherical 
curvature of the wind surface on the photon arrival-time, 
photon frequency Doppler boost, and relativistic beaming. 
Equations (38) and (39) establish the initial distribution of 
the shock-accelerated pairs and the magnetic field, for the 
wind dynamics described in §2. We track numerically the 
evolution of the electron distribution (subject to radiative 
and adiabatic losses) in each shell of external matter, af- 
ter it is accelerated and pair-loaded by the GRB front and 
compressed and heated by the collision with GRB ejecta. 
The radiative and adiabatic losses are also included in the 
calculation of the afterglow dynamics. 



Figure 3. Thick lines: afterglow optical light-curves produced for 
various particle densities n of the circumburst medium, and for a 
fireball initial Lorentz factor F = 100. Other parameters are as for 
Figures 1 and 2. Thin Jincs: light-curves that would be obtained 
if the pair-wind were ignored. Note that the brightening of the 
afterglow emission during the first few minutes by the pair-wind 
increases with the circumburst density. Filled triangles represent 
upper limits obtained by ROTSE, open triangles are for LOTIS 
limits. 

Figures 3 and 4 illustrate the effect of the pair- wind on 
the early optical afterglow for a few values of the circum- 
burst medium density n and fireball initial Lorentz factor T. 
To understand better this effect, we outline first the proper- 
ties of the optical light-curves expected when the radiative 
cooling of the 7p-electrons and the pair-wind are ignored. 
Then the light-curve would peak at time tp when Vp passes 
through the optical band. Equation (41) gives that 

tp^400Ei%et('l,e'J':l,s (42) 

independent of n or F. The optical fiux at tp is the peak 
flux of the synchrotron spectrum, the corresponding optical 

magnitude being 

Ritp) = 14 - 2.5 \og{Ek,s3nl'^e]^%) , (43) 

for redshift 2 = 1. Note that R{tp) increases with the ex- 
ternal density. The light-curves shown in Figures 3 and 4 
with thin lines are consistent with the results in equations 
(42) and (43) for the lowest densities and Lorentz factors 
considered here. As the magnetic field B and typical elec- 
tron Lorentz factor increase with n and F (eqs.[39] and [38]), 
the 7p-electrons cool faster for larger n and F, losing their 
energy on a timescale shorter than the dynamical timescale 
(~ tdec)- As a result, the electron distribution develops a 
tail at energies lower than jp, the peak of the synchrotron 
spectrum moves at a frequency below i/p, and the optical 
light-curve peaks at an earlier time. The peak magnitude 
remains that given in equation (43). 

If the early afterglow spectrum peaks above the opti- 
cal domain (as is expected for values of the micropliysical 
parameters not far from those used in eq.[41]), then the spec- 
trum softening and increase of the number of radiating parti- 
cles produced by the pair wind should brighten the afterglow 
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Figure 4. Same as in Figure 3, but for various cjccta Lorcntz 
factors r and a fixed particle density of external medium n = 
10'*cm~^. Note that for F = 100 the brightening produced by 
the pair-wind is maximal. For F ^ 300 the wind reaches only a 
small fraction of the maximal Lorentz factor before being swept- 
up by the forward shock driven by the GRB ejecta, the effect of 
the wind on the afterglow light-curve becoming negligible. 



Figure 5. Effect of the burst energy, electron, and magnetic field 
parameter on the afterglow optical emission. The canonical set 
of parameters is: n = 10"^ cm^'^, F = 100, Ei^ = Ej = 10^^ erg, 
z = 1, T = 20 s, Ep = niec'-^, a = 1.5, £e = 0.1, es = 10"*, 
p = 2.5. Varying the burst duration, the peak energy of the GRB 
spectrum, or its slope has little effect on the resulting afterglow 
light-curve. 



optical emission until a time of the order of tpair- Figures 3 
and 4 confirm this expectation and also show that the mag- 
nitude of the brightening induced by the pair-wind depends 
on the external density and cjccta initial Lorcntz factor. As 
discussed above, a denser medium and a higher Lorentz fac- 
tor lead to a faster Icpton radiative cooling, which makes the 
optical light-curve peak at an earlier time. For the afterglow 
parameters considered in Figures 3 and 4, external densities 
larger than about n ~ 10^ cm^** lead to light-curve peak 
times earlier than tpair- In this case, the pair-wind brightens 
not only the rising part of the afterglow light-curve but also 
its peak, the effect increasing with n. Increasing the ejecta 
Lorentz factor has a similar effect; however, for a larger F, 
the premature termination of the wind loading reduces the 
brightening produced by the pair-wind. 

Summarizing the results shown in Figures 3 and 4, we 
can say that, without the pair- wind and for a set of repre- 
sentative afterglow parameters, the optical afterglow peak 
should be brighter than J? = 13 for n ^ 10^ cm~^, should 
increase by one magnitude for an ten-fold increase in the 
external density, but this peak should occur earlier than 100 
seconds if n ^ 10'' cm~^. Further, for F ^ 200, the pair- wind 
brightens the optical afterglow peak of the optical afterglow 
by half of a magnitude for n = 10^ cm~^ and by 1 magni- 
tude for n = 10'' cm~^; however, for F ^ 200, the effect of 
the pair-wind is reduced. 

For the canonical afterglows parameters given in Figure 

2 and n ~ 10* cm^"^, the optical afterglow shown in Figures 

3 is brighter than the upper limit f? ^ 13 at t = 10 s set 
by the ROTSE experiment (Akerlof et al. 2000, Kehoe et al. 
2001) for the prompt optical emission of two GRBs; however, 
a lower density or a higher Lorcntz factor (suppressing the 
pair- wind brightening), would explain the non-detections. 
Similarly, the upper limits ranging from ii = 13 to i? = 15, 
obtained by ROTSE at t = 100 s for three bursts, require 



that n ^ 10^ cm~^ and n ^ 10^ cm"'', respectively, the 
pair-wind having a small effect at such times and for such 
densities. Finally, n < 10^ cm~^ is required by the upper 
limits ranging from 7? = 14 to i? = 16 set at t = 10^ s 
by ROTSE for five bursts (see also Rykoff et al. 2002) and 
by R{t > 1200s) > 14 and R{t < 3000s) > 15 obtained by 
LOTIS for GRBs 971227 and 010921, respectively (Williams 
et al. 1999, Park et al. 2002). 

The above comparison between expectations and ob- 
servations for the early afterglow emission suggest that the 
circumburst density does not exceed n = 10^ — 10^ cm~^, un- 
less the microphysical parameters e s and Ee and the kinetic 
energy Ek have lower values than considered for Figures 3 
and 4. The effect of these parameters is illustrated in Figure 
5 for n = 10^ cm~^ and F = 100. A kinetic energy Ek lower 
than 10*^ ergs or magnetic field parameter es below 10~® 
yield an early afterglow emission too dim to have been de- 
tected with the past capabilities of ROTSE or LOTIS, but 
within their reach in the future. For a lower electron energy, 
past detections were more difficult, but not impossible. We 
note that such low values for Ek and the 7-ray output energy 
Ej obtained by Bloom, Frail & Sari (2001) for bursts with 
known redshifts, ranging from 10^'^'^ to lO^**'^ erg, would im- 
ply high burst efficiencies. Furthermore, the above low values 
for Ee and £b are below the values Se > 0.1 and eb > 
determined by Panaitescu & Kumar (2002) from fits to the 
0.5-100 days broadband afterglow emission of several GRB 
afterglows. 

Therefore, barring a substantial increase of the Se and 
Es parameters from 10^ — 10^ s to 10^ s after the burst, 
the non-detections by LOTIS and ROTSE of GRB optical 
counterparts brighter than _R = 13 16 at t < 3 X 10'^ s imply 
that the medium surrounding the GRB source is less dense 
than n = 100 cm~^ up to a distance r ~ 0.01 pc. 

Figures 6 and 7 present the afterglow light-curves ex- 
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Figure 6. Afterglow light-curves for the same parameters as in 
Figures 1 and 2, but for an external medium of density n = Ar~^, 
as expected if the GRB progenitor is a massive star. The legend 
gives the constant A divided by that corresponding to a mass- 
loss rate of 10~* Mq/jt and a stellar wind speed of 10^ km/s. 
Thick lines are for the afterglow light-curve calculations including 
the pair-wind, while thin lines show light-curves obtained if the 
pair-wind were ignored. Note that the pair-wind brightens the 
afterglow emission during the first minute by up to 3 magnitudes, 
its effect disappearing after the second minute, at the peak. 



Figure 7. Afterglow light-curves for the parameters given in Fig- 
ure 6, for the ejccta Lorcntz factors given in legend and a stellar 
wind with A, = 1. For P = 50 the pair- wind reaches initially a 
Lorentz factor 7± larger than that of the GRB cjecta and a gap 
develops between them, the afterglow emission starting when 7-1- 
has fallen below F and the ejecta sweep-up the pair-wind. As for 
Figure 6, the effect of the pair-wind on the afterglow emission 
subsides after a couple of minutes. 



pected when the circumburst medium is a stellar wind of 
density profile n = Ar~^ ejected by the GRB progenitor, 
as expected in the hypcrnova scenario. In contrast with the 
homogeneous medium case considered in Figures 3 and 4, 
the brightening produced by the pair wind is stronger (2-3 
mags), occurs earlier, and has a weak dependence on the nor- 
mahzation constant A (i.e. the wind density), but a stronger 
dependence on the ejecta Lorcntz factor. For the canonical 
burst parameters used here, optical afterglows brighter than 
some of the ROTSE and LOTIS limits are obtained even for 
tenuous stellar winds corresponding to a mass-loss rate to 
wind speed ratio of lO"*^ M0/yr/(lO^ km/s). 



4 CONCLUSIONS 

Assuming that GRBs are produced by internal shocks at 
cm, a small fraction of the gamma-ray photons 
is converted into electron-positron pairs at a distance of less 
than rpair = 10^® cm from the center of explosion. The frac- 
tional energy carried away by this pair-wind is ~ 10~®n, 
where n is the particle density in cm~^ of the unperturbed, 
circumburst medium. The number of pairs created per ex- 
ternal electron is of order a few hundred and the average 
Lorentz factor of the pairs is a few. 

For plausible values of the fireball kinetic energy, elec- 
tron and magnetic field parameters in the shocked circum- 
burst medium, and for an external density n = 100 cm~^, 
the optical afterglow from the forward shock peaks around 
10 minutes from the burst (when the peak of the synchrotron 
spectrum reaches the optical domain) and at magnitude 
i? ~ 14 (Figure 3). For a tenfold increase in the density. 



the peak brightens by one magnitude (eq.[43]) and occurs 

earlier due to the electron cooling. 

By softening the afterglow spectrum and by increas- 
ing the number of radiating particles, the pair-wind bright- 
ens the rise of the optical light-curve. For densities higher 
than n ~ 10* cm~* the peak optical emission occurs at 
radii whore the pair production is significant, leading to 
a brightening of the optical peak emission (Figure 3) if 
the ejecta Lorentz factor is less than ^ 300 (Figure 4). 
A higher Lorentz factor suppresses the pair-wind acceler- 
ation/enrichment and its brightening of the afterglow emis- 
sion until tpair = (1 + z)rpair/{cT'^) ~ 10(r/300)"^ s, but 
leaves unaltered the afterglow emission at later times, past 
the deceleration timescale. 

For conservative values of the afterglow kinetic energy 
and electron/magnetic field parameters, the upper limits set 
by ROTSE and LOTIS on the prompt optical emission from 
several GRBs at 10^ - 10* s, ranging from 13 to 16 mag, re- 
quire that the density of the circumburst medium within 
0.01 pc of the GRB progenitor is less than 100 cm^*. This 
limit is consistent with the typical densities of 0.1 — 30 cm~* 
inferred at r > lO'^'^ cm from fits to the broadband emis- 
sion of several GRB afterglows (e.g. Panaitescu & Kumar 
2002). Moreover, our results show that a wind-like medium 
around the GRB progenitor, corresponding to a mass-loss to 
wind speed ratio above 10~^ M0/yr/(lO* km/s) yields opti- 
cal counterparts brighter than i? = 13 up to almost 10* s. 
Thus, the wind of a Wolf-Rayet star is too dense and incom- 
patible with the ROTSE and LOTIS upper hmits. 

Figures 3, 4, 6, and 7 show that the brightening of 
the pair-wind of the early afterglow emission is much more 
prominent for a wind-like external medium. This feature is 
present within the first 100 seconds of the optical afterglow 
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and provides a test for the type of circumburst environment, 
using future Swift observations. 

For n ^ 10^ cm~^ all photons with energy above the 
threshold for production are transformed into pairs and 
the GRB spectrum is severely modified. In this case the pair- 
wind becomes optically thick to Thomson scattering and the 
burst may not be detectable. However, an optical counter- 
part resulting from the interaction of the ejecta with the 
pair- wind would be very bright [R ~ 10), unless it is atten- 
uated by the dust within the dense medium. Such fiashes, if 
produced, would be missed by ROTSE and LOTIS, in the 
absence of the GRB emission. 
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